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Theory review of exclusive rare radiative decays 
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I briefly review the theory status of exclusive rare radiative decays. 
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1. Introduction 

Experimental measurements of B — > V7 decays, 
with V a light vector meson such as K* , p, u>, (f>, have 
continued to improve and will become more precise at 
the end of the B-factories and at LHC-b. Given their 
rich CKM phenomenology and potential to constrain 
new physics models, having reliable theory predictions 
for these decays is increasingly relevant. In this talk I 
briefly review the status of this area. 

The starting point is the effective weak Hamilto- 
nian, which for B — ► decays is [l|: 



CiQi + C2Q2 + CiQi 



i=3 



= % E W 

* p—u,c 

(1) 

where = V* q V p b- The operators with the largest 
Wilson coefficients are the four-quark operators Q\ 
and Q2, which read 



Qi = (qp)v-A{pb)v-A 
Q P 2 = {qiPj)v-A(Pjbi)v-A, 



(2) 



and the electromagnetic and chromomagnetic penguin 
operators Qi and Qg, which are 

Qi = - e -^-qa^[l + 75 ]bF, u , 



Q s = - 



8tt 
gm b (p) 
8tt 2 



qa^[l + l5 ]T a bG%. (3) 



Here q — - d or s, and the convention for the sign of 
the couplings corresponds to the covariant derivative 
iDp = id^ + eQfA^ +gT a A a ^ with and A^ repre- 
senting the photon and gluon fields respectively, and 
Q e = — 1 etc. The factor mb(p) is the MS mass of the 
b quark. 

The main theoretical challenge is to evaluate the 
hadronic matrix elements of the operators in the effec- 
tive weak Hamiltonian. Common ways of doing this 
include the QCD factorization @, H, 0| and pQCD 
[E [H approaches. Both of these rely on the fact that 
ct s (jnb) and Aqcd /w-b can be considered as small ex- 
pansion parameters. In this talk I will focus exclu- 
sively on QCD factorization methods, including in this 
category also strategies which supply additional infor- 
mation on 1/mb-suppressed contributions from QCD 
sum rules 0, , or apply a form of renormalization- 
group (RG) improved perturbation theory based on 
soft-collinear effective theory (SCET) [UEJ. 



The remainder of the talk is organized as follows. 
Section[2]introduces the QCD factorization formalism, 
Section [3] gives two sample applications, and Section 
0] reviews recent results in higher-order perturbative 
corrections to the hard-scattering kernels. This is fol- 
lowed by a short discussion of 1/ mj, power corrections 
and endpoint divergences in Section [5l using isospin 
violation in B — > K*^ decays as an example, and a 
summary in Section O 



2. QCD factorization 

QCD factorization is the statement that in the 
heavy-quark limit the hadronic matrix element of each 
operator in the effective weak Hamiltonian can be 
written in the form 



(Vl\Q t \B) 



tl+ 



(4) 



O 



A 



QCD 



where the ® stand for convolution integrals. (Some- 
times this same formula is written in a different but 
equivalent form, using the tensor QCD form factor 
Tf^ v instead of the SCET soft function Cvi-) The 
soft function (v ± , the meson decay constants fv ± , F, 
and the light-cone distribution amplitudes (LCDAs) 
</>5, are non-perturbative but universal objects. At 
present, numerical values for these objects are taken 
from light-cone sum rules. The hard-scattering ker- 
nels i*' 11 can be calculated as a perturbative series in 
a s . The t] are referred to as "vertex corrections", and 
the t] 1 are referred to as "spectator corrections" . Only 
the hard-scattering t\ is non- vanishing at zeroth order 
in a s . Corrections from this and the other operators 
appear at order a,, and have been known completely 
for some time @, y, 0] • I shall refer to these a s correc- 
tions as next-to- leading order (NLO), even though the 
hard-spectator kernels t\ l first start at this order. In 
Sectional describe some recent results for the NNLO 
kernels. 

An all orders proof of the QCD factorization for- 
mula (|U) was performed in Q, using the technology 
of SCET. In the effective-theory approach the hard- 
scattering kernels are short distance Wilson coeffi- 
cients of operators whose hadronic matrix elements 
define the (y ± and the LCDAs. For the vertex term 
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this makes little practical difference, but for the hard- 
spectator term one can show that the hard-scattering 
kernels t] 1 can be further factorized into the form 

t] l (u,u),fii) = [C? 1 (fj, h )®U(fi h> iJ li )]®j ± (Lii)- (5) 

The hard coefficients Cf 1 contain physics at the scale 
TOf,, while the jet function j± is independent of the op- 
erator and contains physics at the intermediate scale 
U mf,AQCD ~ 1.5 GcV. The evolution factor U is de- 
rived by solving the RG-equations in the effective the- 
ory. Because these equations are non-local, the evolu- 
tion factor appears in a convolution with the Wilson 
coefficient Cf 1 . In the limit where m\ ^> m^ApcD, 
this evolution matrix resums perturbative logs in the 
ratio Aqcd/™&, and allows one to evaluate the hard 
coefficients and the jet function at their natural scales 
/i/i ~ mb and /i, ~ mbKqpD- If one does not wish to 
do this, then making the choice /i/j = fa sets the evo- 
lution factor U to unity, and one recovers the original 
QCD factorization formula ([4]). So for these decays, 
the "SCET approach" is just QCD factorization, plus 
the added opportunity to perform a further scale sep- 
aration in the hard scattering kernels t] 1 . 

It is necessary to keep in mind that the factor- 
ization formula Q is valid only up to 0(Aqcd/?ti6) 
power corrections. Since SCET is an effective the- 
ory which sets up a systematic expansion in a s and 
Aqcd/^6, it has the potential to provide a complete 
classification of the subleading terms in 1/mj,. How- 
ever, there has been no serious attempt to do this, 
because the generalization of the leading-order fac- 
torization formula to corrections in 1 /m& presumably 
contains many more non-perturbative objects. Even 
more troublesome is the fact that, in the cases where 
power corrections been calculated, the convolution in- 
tegrals over momentum fractions do not always con- 
verge. These "endpoint divergences" are at present a 
principle limitation on the entire formalism. 

In the absence of a comprehensive theoretical frame- 
work to deal with the 1 / suppressed corrections, it 
has become common practice to focus on the contri- 
butions which are believed to be large, or which play 
an especially important role for phenomenology. In 
this talk I will discuss the power corrections which 
are important for the calculation of isospin violation 
in B — > K*j decays as performed in ; ll|. This is 
sufficient to describe the problem of endpoint diver- 
gences. Other power corrections, due to long-distance 
photon emission and soft-gluon emission from quark 
loops, were calculated within the context of light-cone 
sum rules in Q . 

3. Sample applications 

The formalism above can be used to calculate many 
different observables, see for instance 

13111113. In 



this section I briefly touch on two examples, taking 
results from the recent studies 0, @| for simplicity. 
The first is the determination of the ratio of CKM 
matrix elements |Vjd/Vt s | from the ratio of branching 
fractions in B — > p7 and B — > K*^y decays. This is 
an independent check on the measurements from B s 
oscillations The most recent theory results are 



~ B(B^K*j) \V ts \* W 
x (0.75 ± 0.11(£) ± 0.02(UT param., 0(l/m b ))) , 

which combined with February 2007 HFAG data leads 
to |Vtrf|/]Vta| = 0.192 ± 0.014(th) ±0.016(exp). This 
central value is compatible with that from B s oscil- 
lations, although the experimental errors are signif- 
icantly larger and the theory errors about twice as 
large. The dominant theory error is in the form factor 
ratio £ = T®^ K " /T^ p ', and according to the recent 
review [l6| this should go down by about a factor of 
two with improved lattice results for the tensor decay 
constants. It should be mentioned that the experi- 
mental measurements of the ratio in ^ combine data 
from p and u> decays under the assumption of exact 
isospin asymmetry, while the theory calculations can 
actually predict the magnitude of isospin breaking. 
Im pro ved experimental measurements such as those 
in [l7l [l8j will make this isospin averaging unneces- 
sary. 

As a second application I quote results for the 
isospin asymmetry in B — > K*^/ decays. This asym- 
metry was calculated in the context of QCD factor- 
ization in [ll| and the numbers updated in [8|, which 
gives the result 

A {K *. = r(B° ^K*°j)-r(B- ->ir- 7 ) 

/l ' T{B° -> K*°-y) +F(B- -> K*~j) 

= (5.4 ±1.4)%. (7) 

The current result by HFAG in [l9[ is (3±4)%. Results 
from [Tlj showed that this asymmetry is particularly 
sensitive to the penguin operator Qq. Improved ex- 
perimental measurements of this asymmetry can thus 
provide constraints on new physics. The theory pre- 
diction above, however, does not include any a s cor- 
rections to the perturbative hard-scatterin g ke rnels, 
even though some of these were calculated in [ll| . The 
reason for this is the presence of endpoint divergences, 
a topic discussed in more detail in Section [5] 



4. NNLO perturbative corrections 

Recently, a set of the NNLO corrections to the ver- 
tex and spectator kernels were obtained in [lOj. In 
this section I summarize the new results from that 
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work and also what remains to be done to complete 
the NNLO calculation. 

Consider first the vertex corrections and the hard- 
scattering kernel t 1 . This hard-scattering kernel is 
most easily obtained as a matching coefficient in 
SCET; details were given in [|, The only tech- 
nical point I mention here is that the matching co- 
efficients are independent of the external states used 
in the calculation, and it is possible to use the par- 
tonic matrix element (qj\Qi\b) in the matching. The 
SCET calculation is trivial for on-shcll quarks, be- 
cause the loop corrections are given by scaleless in- 
tegrals which vanish in dimensional regularization, so 
the main challenge is to calculate the partonic matrix 
elements in full QCD. However, these are just vir- 
tual corrections to the b — > (77 process needed also for 
inclusive B — > X s j decay, and are currently known 
very accurately due to the efforts of many people 
(see [l(J H3] for a complete list of references). Us- 
ing those multi-loop calculations, t^ s can be obtained 
completely to NNLO. For the four-quark operators 
Qi,2, on the other hand, results to NNLO are known 
only in the large- (3q limit, obtained by calculating the 
a 2 s Cpnf terms and replacing nf — > — 3/2/?o, where 
f3a = 11 — 2/3n./ in QCD. The remaining terms re- 
quire to calculate a large number of three-loop graphs 
depending on the ratio m 2 /m 2 and containing imagi- 
nary parts. Given that this is also a missing piece in 
the B — ► X s j calculation, it is likely that this will be 
done in the near future. 

I now turn to the spectator corrections and the hard 
scattering kernels t\ . Again in this case their calcu- 
lation is most conveniently formulated as a matching 
calculation in SCET. In contrast to the vertex correc- 
tions, however, the Feynman diagrams needed in the 
matching calculation have no analog in inclusive de- 
cay, where interactions with the spectator quarks do 
not contribute at leading power in l/m&. Therefore, 
results for the t\ l cannot simply be extracted from ex- 
isting calculations and must be obtained from scratch. 
The one exception is the a 2 contribution to t l 7 l , which 
can be taken from studies of heavy-to-light form fac- 
tors carried out in l2ll |22| . The a 2 contribution to ig 1 
was calculated in |l0|. The most important missing 
pieces are from <3i,2, but these are also the hardest to 
calculate, since they involve two-loop graphs depend- 
ing on m 2 /m 2 in addition to the momentum fraction 
u of the quark in the V-meson. 

A natural place to explore the numerical impact 
of the NNLO corrections is on branching fractions in 
B — > K*-f decays, since for these decays the annihi- 
lation topology is CKM suppressed. The branching 
fraction is 



-in 



(8) 

where the A v (Ahs) are the amplitudes for the ver- 



tex (hard spectator) corrections. These are separately 
RG-invariant and so can be studied individually. For 
the vertex corrections, the ratio of the NNLO ampli- 
tude to LO amplitude for the default set of parameters 
in pi] is 



A NNLO 



= 1 + (0.096 + 0.057z) [a s ] 
+ (-0.007 + 0.030i) [a 2 ] , 



where the first term in parentheses is the NLO (a s ) 
correction and the second term the NNLO (a 2 ) cor- 
rection. One sees that the real part of the NNLO cor- 
rection is extremely small. It is instructive to further 
split the above amplitude into the pieces originating 
from the various operators Qi . In that case the result 
reads 



NNLO 



1 = 



(9) 



(.26 + .03i) [Qi] - .18 [Q 7 ] + (.02 + .02i) [Q 8 ] [a 



+ ( (.07 + .02i) [Qi] - .08 [Qr] + (-002 + Mi) [Q 8 



From this one sees that the NNLO contributions 
are so small because of a large cancellation between 
the Qi and Q7 contributions. It is an open question 
whether this cancellation will persist when results for 
Qi are obtained beyond the large-/3o approximation 
used above. This will be discussed in more detail be- 
low, where results for branching fractions are given. 
For the hard spectator amplitude, one finds 

/(NNLO 

hs LO = (0.11 + 0.05i) [a s ] + (0.03 -I- O.Olz) [a 2 s ]. 

(10) 

Unlike the case of the vertex corrections, the indi- 
vidual contributions from the different operators are 
rather small at NLO and especially NNLO. The exact 
numbers, including a split into components at the jet 
and hard scales in ([5]), can be found in [Tol | . 

Putting the amplitudes together, one can find the 
branching fractions and their uncertainties. Includ- 
ing isospin and SU(3) breaking from meson masses, 
lifetimes, and U-meson distribution amplitudes, (ic| 
estimated 



B(B+ 
B(B° 



B(B S -► 7 ) 



(4.6 ± 1.4) x 10~ 5 , 
(4.3 ± 1.4) x 10~ 5 , 
(4.3 ± 1.4) x 10~ 5 . 



(11) 



The magnitude of various sources of uncertainty can 
be found in [10]. Most significant is about a 25% 
uncertainty due to the soft functions Cvi 1 which is 
expected. Somewhat surprising is that the a 2 cor- 
rections from Qi can make a large impact, and since 
they are taken only in the large-/?o limit the true re- 
sult is still rather uncertain. For instance, assigning a 
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Figure 1: Feynman diagrams contributing to isospin vio- 
lation in B — > 7^*7 decays. 

100% uncertainty to the a 2 s piece from Q\ in ([9]), one 
finds an uncertainty of about ±0.5 in the branching 
fractions. Also, because the CpUf terms do not fix 
the perturbative definition of the charm-quark mass 
in the a s contribution in ((9|), one can formally use 
values ranging from a low MS mass to a high pole 
mass, which can make a large numerical difference, 
about ±0.4 in the range used in [l(|- The message to 
be gained from this is that the full NNLO corrections 
from the set Q\ t i also need to be calculated, in order 
to get the reduced perturbative uncertainty expected 
from a higher-order calculation. 

5. Isospin violation and endpoint 
divergences 

Calculations of interesting observables such as 
isospin violation in B — ► K*"f decays require to in- 
clude l/nib suppressed contributions to the factoriza- 
tion formula. Particularly important corrections come 
from topologies such as those in FigureHJ which I have 
taken from [ll( . The square vertex indicates an inser- 
tion of Qi . . . Qe for the graph on the left, an insertion 
of Qs for the graph in the center, and an insertion of 
Qi.2 for the graph on the right. The crosses denote al- 
ternate attachments of the photon line. Note that the 
center and right graphs are a s corrections compared 
to the graphs on the left. The contributions from the 
graphs in the figure, as well as some a s terms linked 
to the graph on the left by RG-invariance, were cal- 
culated in [ll[ . To finish the calculation to this order 
would require to calculate the a s corrections to the 
graphs on the left. This calculation is in progress [23j . 
However, once completed, the calculation will still face 
the problem that a straightforward application of the 
QCD factorization formalism leads to endpoint diver- 
gences in the convolution integrals. For example, the 
center graph in Figure [1] was shown in (ll| to contain 
the integral 

x ± = f\tfM 1 *^ s , (12) 

which is equal to infinity under the conventional as- 
sumption that the LCDA vanishes as 1 — u in the 



endpoint. The quantity X±_ is multiplied by the Wil- 
son coefficient Cs and was estimated by introducing 
an IR cutoff in [ll| and found to be small, but clearly 
this is a conceptual problem. Moreover, there is no 
guarantee that similar endpoint divergences will not 
appear for the other operators in the weak Hamilto- 
nian, once the full set of one-loop corrections to the 
graphs on the left in Figure[T]are included. Therefore, 
it is fair to say that a better understanding of power 
corrections is needed, in order to achieve a consistent 
framework for the calculation of quantities such as 
isospin violation. 



6. Summary 

Rare radiative B — > V~f decays are of increasing 
interest as experimental measurements become more 
precise. QCD factorization and SCET have provided 
a theoretical framework which can be used to calcu- 
late observables for these decays to leading order in 
1/mb- The perturbative hard-scattering kernels have 
been known at NLO in a s for some time, and recently 
a set of NNLO results have been obtained for the 
operators Qi, Qt, and Q$, although results for the 
four-quark operators Q\ and Q2 are not yet complete. 
The interesting observables such as isospin asymme- 
tries and branching fractions in the b — > dj modes are 
sensitive to power corrections in l/rrife. Some of these 
have been estimated in the framework of light-cone 
sum rules, but a systematic treatment in SCET or 
QCD factorization is still missing due to the presence 
of endpoint divergences; a solution to this problem 
would be a much desired advance in this field. 
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